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Abstract 

We present algebraic construction of Darboux matrices for 1+1- 
dimensional integrable systems of nonlinear partial differential equa- 
tions with a special stress on the nonisospectral case. We discuss 
different approaches to the Darboux-Backlund transformation, based 
on different A-dependencies of the Darboux matrix: polynomial, sum 
of partial fractions, or the transfer matrix form. We derive symmetric 
A^-soliton formulas in the general case. The matrix spectral parameter 
and dressing actions in loop groups are also discussed. We describe 
reductions to twisted loop groups, unitary reductions, the matrix Lax 
pair for the KdV equation and reductions of chiral models (harmonic 
maps) to SU (n) and to Grassmann spaces. We show that in the 
KdV case the nilpotent Darboux matrix generates the binary Dar- 
boux transformation. The paper is intended as a review of known 
results (usually presented in a novel context) but some new results 
are included as well, e.g., general compact formulas for A^-soliton sur- 
faces and linear and bilinear constraints on the nonisospectral Lax 
pair matrices which are preserved by Darboux transformations. 
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1 Introduction 



A 1+1- dimensional integrable system can be considered as integrability con- 
ditions for a linear problem (a system of linear partial differential equations 
defined by two matrices containing the spectral parameter) , see for instance 
|56j . The Darboux-Backlund transform is a gauge- like transformation (de- 
fined by the Darhoux matrix) which preserves the form of the linear problem 
[H [221 123 Sni [66] . All approaches to the construction of Darboux matrices 
originate in the dressing method [561 [SHI HH [82] . 

The paper is intended as a presentation of Darboux-Backlund transforma- 
tions from a unified perspective, first presented in [131 E]. The construction 
of the Darboux matrix is divided into two stages. First, we uniquely char- 
acterize the considered linear problem in terms of algebraic constraints (the 
divisor of poles, loop group reductions and other algebraic properties, e.g., 
linear and biblinear constraints). Then, we construct the Darboux matrix 
preserving all these constraints. Using general theorems, including those 
from the present paper, one may construct the Darboux matrix in a way 
which is almost algorithmic. 

The paper is intended as a review of known results but some new results 
are also included. We discuss in detail elementary Darboux transformation 
(Darboux matrix which has a single simple zero), symmetric formulas for 
Darboux matrices and soliton surfaces (in the general case), and loop group 
reductions for polynomial Darboux matrices. Two examples are discussed in 
detail: the Korteweg-de Vries equation and chiral models (harmonic maps). 

The part which seems to be most original contains the description of 
linear and bilinear invariants of Darboux transformations. We prove that 
multilinear constraints introduced in [H] are invariant with respect to the 
polynomial Darboux transformation (also in the nonisospectral case). Taking 
them into account we can avoid some cumbersome calculations, our construc- 
tion assumes a more elegant form and, last but not least, we do not need any 
assumptions concerning boundary conditions. 

Another important aim of this paper is to show similarities and even an 
equivalence between different algebraic approaches to the construction of the 
Darboux matrix. This is a novelty in itself because sometimes it is difficult 
to notice connections between different methods. The existing monographs, 
even the recent ones, focus on a chosen single approach, compare [211 [23 [13 
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We consider a nonlinear system of partial differential equations which is 
equivalent to the compatibility conditions 

U^,,-U,,^+[U^,U,] = Q , (l^/i<z/^m), (1.1) 

for the following system of linear equations (known as the Lax pair, at least 
in the case of two independent variables) 

^,,= f/,^, (z/ = l,...,m) , (1.2) 

where n x n matrices U,y depend on a;^, . . . , x*" and on the so called spectral 
parameter A (and, as usual, \E',,y = /dx^ , etc.). We assume that is also 
a matrix (the fundamental solution of the linear system (11. 2p ). We fix our 
attention on the case m = 2 (although most results hold for any m) and 
shortly denote by x the set of all variables, i.e., x = [x^, . . . , x™). 

1.1 Isospectral and nonisospectral Lax pairs 

Let us recall that the most important characteristic of the matrices Ui, U2 
is their dependence on the spectral parameter A. In the typical case Uiy are 
rational with respect to A. Actually we will consider a more general situation. 
We assume that the Lax pair is rational with respect to A, and 

• "isospectral case": A is a constant parameter, 

• "non-isospectral" case: 

X,^=L^{x,\), (z/ = 1, . . . ,m) , (1.3) 

where Liy are given functions, rational with respect to A (this case 
reduces to the isospectral one for Li^{x, A) = 0). 

Remark 1.1. The differential equations U.3\) are of the first order, so their 
solution A = A(a;, () depends on a constant of integration ( which plays the 
role of the constant spectral parameter. 

The solution of the system (11.31) exists provided that compatibility con- 
ditions hold, for more details see [H]. In general A = A{x,() is an implicit 
function, although in many special cases explicit expression for A can be 
found, compare [UlIIllES])- 
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1.2 The Darboux-Backlund transformation 



The application of the dressing method to generate new solutions of nonlinear 
equations "coded" in (11.11) consists in the following (see [561 [79l [82] ) . Suppose 
that we are able to construct a gauge-like transformation \l/ = (where 
D = D{x,X) will be called the Darboux matrix) such that the structure of 
matrices U^, 

U^ = D,,D-' + DU,D-' , (z/= l,...,m) , (1.4) 

is identical with the structure of the matrices U^. The soliton fields entering 
f/jy are replaced by some new fields which, obviously, have to satisfy the 
nonlinear system (11.11) as well. 

Remark 1.2. The Darboux transformation should preserve divisors of poles 
(i.e., poles and their multiplicities) of matrices U^. This is the most impor- 
tant structural property of to he preserved. The second important property 
is the so called reduction group, see Section\^ 

For any pair of solutions of (11.11) one can "compute" D := The 
crucial point is, however, to express D solely by the wave function \E' because 
only then one can use D to construct new solutions. Such D is known as 
the Darboux matrix [lOl HSl HB]. The Darboux matrix defines an explicit 
map 5* 1-^ 5*, where 5" is the set of solutions of the linear problem (11.21) . The 
construction of the Darboux matrix is based on the important observation: 

Remark 1.3. The Darboux matrix can be expressed in an algebraic way by 
the original wave function . 

By the "original wave function" we mean one before the transformation. 
In fact, it is rather difficult to find special solutions of the linear problem. 
Usually very limited number of cases is available. However, knowing any 
solution \l/ = \l/(a;. A) and the Darboux matrix one can generate a sequence 
of explicit solutions. Starting from the trivial background (x-independent 
and mutually commuting Uu) we usually get the so called soliton solutions. 

1.3 Equivalent Darboux matrices 

It is quite natural to consider as equivalent Darboux matrices which pro- 
duce exactly the same transformation (II. 4p of matrices Ui, of a given linear 
problem. 
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Remark 1.4. The linear problem U.2\) is invariant under transformations 
\& \I/Co (for any constant nondegenerate matrix Cq = Co{X)). 

Therefore Darboux matrices D and D' are equivalent if there exists a 
matrix C such that = D'^^C (for any ^). Thus C should commute with 
\E' what, in practice, means that C = /(A) G C. 

Remark 1.5. The matrix D' = f{X)D, where f is a complex function of X 
only, is equivalent to D. 



1.4 Soliton surfaces approach 

Given a solution \E' = \E'(x, A), where A depends on x and (, we define a new 
object F by the so called Sym-Tafel (or Sym) formula: 

F = ^-'^,^, (1.5) 

If \E' assumes values in a matrix Lie group G, than (for any fixed () F describes 
an immersion (a "soliton surface" ) into the corresponding Lie algebra [7T1 [71] . 
Soliton surfaces are a natural frame to unify a variety of different physical 
models like soliton fields, strings, vortices, chiral models and spin models [72] . 
In the framework of the soliton surfaces approach one can reconstruct many 
integrable cases known from the classical differential geometry [71 |15l [161 [H] ■ 
The Darboux-Backlund transformation for soliton surfaces reads 

F = F + ^-^D-^D,t;^ (1.6) 

where Z),^ = A,^ D,\. The equivalent Darboux matrices yield the same soliton 
surfaces. Indeed, if we take D' = fD, then 

F = F + ^^-^D'-^D',^ ^^J = F + ^ + ^-^D-^D,^ ^ , (1.7) 

i.e., surfaces corresponding to D and D' differ by the constant (In/),^. 

In order to illustrate usefulness of the geometric approach we present the 
following theorem |18j . 

Theorem 1.6. We assume that Ui,U2 are linear combinations ofl, X and 
X~^, with x-dependent su{2)-valued coefficients, and Uv{—X) = EoUu{X)Eq^ 
(where Eq e su{2) is a constant matrix). Then F given by the Sym formula 
M.5\) is (in the isospectral case) a pseudo spherical (i.e., of negative Gaussian 
curvature) surface immersed in su{2) ~ M'^. In the nonisospectral case the 
same assumptions yield the so called Bianchi surfaces. 
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We point out that surprisingly few assumptions (restrictions) on the spec- 
tral problem leads to the very important class of pseudospherical surfaces. 
It is easy to assure the preservation of these restrictions by the Darboux 
transformation. 

Darboux transformations usually preserve many other constraints (e.g., 
linear and bilinear invariants discussed in Section [S]) what leads to the preser- 
vation of some geometric characteristics (e.g, curvature lines) and to a specific 
choice of coordinates and other auxiliary parameters. 

2 Binary Darboux matrix 

In this paper by the binary Darboux matrix we mean one pole matrix with 
non-degenerate normalization 

D = Af(^I + ^^l^P^ ' P^ = P^ detAT^O, (2.1) 
such that its inverse has the same form: 

D-'=(l + l^l^p] M-^ . (2.2) 

V A - /ii y 

Here Ai,/ii are complex parameters (which can depend on x in the non- 
isospectral case), P = P{x) is a projector matrix (P^ = P), and A/" = //{x) 
is the so called normalization matrix. 

2.1 Binary or elementary? 

The name "binary" for Darboux matrices of the form (12.11) is rather tentative, 
because binary Darboux transformations were introduced in another context 
(compare P6| l83^ j ). The "classical" binary transformation corresponds to the 
degenerate case of (12.11) when /ii —> Ai (see Section |44|) . i.e., 

where = (the so called nilpotent case, see [E]). Therefore, we use this 
notion in an extended sense. However, it seems to be compatible with under- 
standing binary Darboux transformation as a composition of an elementary 
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Darboux transformation and a Darboux transformation of the adjoint linear 
problem [IBl \57l EZj- In the case of the Zakharov-Shabat spectral problems 
(11.21) the adjoint spectral problem is given by 

<^U,, (2.4) 

and one can easily check that $ = solves the adjoint spectral problem. 
The general solution of (12.41) is $ = \E'~^C, where C is a constant (i.e., x- 
independent) matrix. We will see in Section 15.11 that the binary Darboux 
matrix can be expressed by a pair of solutions: one solves the spectral prob- 
lem f ll.2p . and the second one solves the adjoint problem (12.41) . 
The matrix (12.11) is equivalent to the linear in A matrix D 

^ = Ar(A- Ai + (Ai-/ii)P) . (2.5) 

Darboux matrices linear in A are sometimes referred to as "elementary", 
see [60]. Indeed, iterating such transformations we can get any Darboux 
transformation with nondegenerate normalization. However, we reserve the 
name "elementary" for matrices which are not only linear in A but have a 
single zero (see Section H]), or even a single simple zero. The polynomial form 
(12.51) of the binary Darboux matrix has two zeros: Xi,fii. The sum of their 
multiplicities is n. Therefore these zeros are simple only in the case n = 2. 



2.2 Sufficient conditions for the projector 

Assuming that f/j, are regular (holomorphic) at A = Ai and A = /ii, and 
demanding that Uu (expressed by (II. 4p ) have no poles at A = Ai and X = fii 
as well, we get the following conditions (for vanishing the corresponding 
residues), compare [H] : 

Po(-9, + f/,(Ai))o(J-P) = 0, 

(2.6) 

{I-P)o{-d, + UM)oP = 0, 

Ai,^= L^(x, Ai) , fii,^= L^{x,ni) , (2.7) 

where the circles mean composition of linear operators and L,^ are defined by 
(11.31) . Note that for any operators A, B we have: 

A o P = ^ imP C ker A . (2.8) 

Indeed, {A o B)ip = for any vector i.e. A{Bip) = what means exactly 
that Bip G ker A. On the other hand, any element of imB is of the form B(p. 
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Remark 2.1. The assumption that are regular at X = Xi and X = jJi 
(assumed throughout this paper) is essential. Relaxing this requirement we 
can get solutions different from those obtained by the standard Darboux-Bdck- 
lund transformation. Solutions of this kind ( unitons ) have been found in the 
case of harmonic maps into Lie groups fn\j , see also JWTl 

If the system (11.31) has the general solution A = A(a;, (), then the equations 
(12.71) can be solved in terms of the function A: 

Ai=A(x,Ci), fii = A{x,C[), (2.9) 

where Ci;Ci ^I's constant parameters, compare [HEl]. Taking into account 
kerP = im(J — -P), we easily show that the system (12. 6p is equivalent to: 

(-9^ + [/^(Ai))kerP C kerP , 

(2.10) 

{-d, + U,{fii))imP C imP . 

Now we easily see that the conditions (12.101) are satisfied by the projector 
defined by the Zakharov-Shabat formulas (compare [551 [T^f 152]): 

kerP = ^iX,)Vker , imP = ^(/xi)l^.^ , (2.11) 

where \E'(Ai) = \l'(x, Ai), ^^(/Ui) = \E'(x, /ii) and Vker and Vim are constant 
vector spaces such that V^er © Vim = V. Indeed, in this case, by virtue of 
(11.21) . the left-hand sides of (I2.10p are simply equal to zero. 

Taking into account that any projector P can be expressed explicitly by 
its kernel and image, P = (imP, 0)(imP, ker P)~^, we can summarize the 
above discussion as follows. 

Proposition 2.2. The transformation ( [i.^| ) with D given by Ii2. where 

p = {^{^ll)Vm, ^){^{^il)Vim. ^(Ai)\4er)-' , (2.12) 

preserves the divisors of poles of matrices Up. 

The formula (12.121) yields a sufficient condition for P to generate the 
Darboux matrix. It is interesting to find also necessary conditions. Therefore, 
we will try to obtain the most general solution of (I2.10p . 
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2.3 The general form of the binary Darboux matrix 

It is convenient to represent vector spaces in a matrix form. Namely, if 
wi, . . . ,Wk span a vector space V, then we can identify V with the matrix 

V = {wi,...,Wk) . (2.13) 

This matrix has k columns (wi, . . . , w^^) and n rows {n = dim V"). 

Note that because of the freedom in choosing a basis in the vector space 
there are many matrices representing the same vector space. If a^- are coef- 
ficients of a. kxk non-degenerate matrix A, then the vectors 

k 

1=1 

form another basis in V which can be represented by the matrix 

V'^{w[,...,w',)^VA. 

The matrices V and V (for any non-degenerate A) represent the same vector 
space and, in this context, are considered as equivalent ones. 

The space of /c-dimensional subspaces of an n-dimensional vector space 
over C is known as Grassmannian Gk,ni'C) (of course, considering real vector 
spaces we have real Grassmannian Gk^ni^))- The elements of the Grassman- 
nian are classes of equivalence of kxn matrices with respect to the equivalence 
relation: y ~ y if there exists kxk matrix A (det A^ 0) such that V — VA. 

Therefore using the same notation for the vector space and the matrix 
representing it, one should remember about this equivalence. In particular, 
in order to show that some vector spaces W and V are identical, one has 
to consider the equation W = VA with an arbitrary non-degenerate A. In 
the similar way one can check whether is a subspace of V (of course a 
necessary requirement is dim V > dim W) . 

Proposition 2.3. Let V,W be vector spaces, k' = dimW^ < dimV = k. 
Then W dV if and only if there exists kxk matrix B such that W = VB. 

Proof: liW <ZV, then there exists a basis of V such that its first k' vectors 
span W . We represent V by vectors of this basis, i.e.. we choose A such that 
w^, . . . , w'l^i span W . Finally, we put B = A diag(l, . . . , 1, 0, . . . , 0). □ 
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Note that formally W and V belong (in general) to different Grassman- 
nians. But if det B = 0, then the columns of VB are linearly dependent and 
VB can be treated as an element of a Grassmannian of lower dimension. 

We proceed to solving the system fl2.10p . Assuming det\l/(Ai) 7^ and 
det^'(/ii) 7^ 0, we can always put kerP = 'if{Xi)V, imP = where 
W, V are some vector spaces (in general x-dependent). Substituting to fl2.10l) 
we have: *(Ai)l^,^ C ^(Ai)V^, ^(/ii)iy,^ C ^(/ii)^^. Hence, 

By Proposition [231 "we rewrite V,i, C as V,^, = VB^ for some B^, (which have 
to satisfy appropriate compatibility conditions), and analogical equations for 
W. Taking into account the freedom of changing the basis when changing x: 
V = VA (det A^O),we obtain: 

{V'A-'),,= V'A''B, 

Therefore, choosing A such that 

{A-'),,= A-'B, , (2.14) 

we obtain V',^, = 0, i.e. there exists an x-independent basis in V (the same 
conclusion holds for W). The solution of ( I2.14p exists because B^ satisfy 
the compatibility conditions mentioned above. Thus we have shown that the 
formulas (12. lip give the most general solution of fl2.10p . 



3 Polynomial Darboux matrices: general case 

In this paper we consider only rational Darboux matrices [n x n matrices 
with coefficients which are rational functions of A). 

Remark 3.1. In the isospectral case, every rational Darboux matrix is equiv- 
alent to a polynomial Darboux matrix 

N 

D = J2Tk{x)X''-' . (3.1) 

Indeed, it is enough to multiply given D by the least common multiple of all 
denominators. The obtained polynomial will be denoted by D. 
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Another equivalent form of D is a polynomial in A~^, obtained from D{X) 
by dividing it by A^. In some cases this polynomial is more convenient that 
D because it is analytic at X = oo. 

In the nonisospectral case the least common multiple of all denominators 
depends on x. Therefore, any rational Darboux matrix is equivalent to some 
polynomial matrix up to a scalar x-dependent factor. 



3.1 The determinant of the Darboux matrix 

The trace of a quadratic matrix is defined as the sum of diagonal elements of 
this matrix. Both the trace and the determinant are invariant with respect 
to similarity transformations: Tr{BAB~^) = TiA, det{BAB~^) = det A. 

Theorem 3.2 (Liouville). If'^-,v = U^"^ , where v is fixed and Uy = Uy{x) is 
given, then 

(det^),^= Trt/^det^ . (3.2) 

This theorem is well known as the Liouville theorem on Wronskians, see, for 
instance, [2]. 

Applying the Liouville theorem to the Darboux transform \E' = D'^ we 
get (det D det = Tr[/i^ det D det Hence, using once more (13.21) . we 
obtain: 

(detD),, 



deiD 



Tr(?7,) - Tr(?7,) . (3.3) 



Remark 3.3. We usually consider traceless linear problems (Ti 11,^ = for 
V = l,...,m). In such case det D has to be constant (i.e., det D does not 
depend on x). Therefore, in the isospectral (and traceless) case detD can 
depend only on X and all its zeros are constants. 

In the nonisospectral case the situation is more complicated because A 
depends on x. However, it is still possible to obtain a strong general result 
characterizing zeros of detD. 

Theorem 3.4. We consider a polynomial Darboux matrix D for a non- 
isospectral linear problem U.^) with X satisfying U.3\) . If detD{Xk) = 
and matrices Uu are regular at Xk, then 

Xk,u= Ly{x,Xk) , (3.4) 

i.e., Xk = A{x,(k), where (k = const. 
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Proof: The determinant of the polynomial -D(A) has a finite number of roots 
(x-dependent, in general). We denote them by A^, k = 1, . . . ,K, and their 
multiplicities by rrik- Note that mi + m2 + . . -mK = nN, where is the 
degree of the polynomial -D(A) and n is the order of the matrix D. Thus 



K 



detZ)(A) =/iJJ(A-Afc)™'= , (3.5) 



fc=i 



where h = h{x) and Afc = Afc(x). Taking into account (II. 3p we compute 



K 

(detD),^ h,^ , / -^i.(a;, A) - Afc, 



+ 1^ U. — . (3.6) 



det D h ^ \ A — Afc 

fc=i 

The equation (13. 3p with regular at A^ implies that the right-hand side of 
(13.61) should have no poles. Therefore residua of (13.61) at A = A^ vanish what 
implies (13. 4p . The x-dependence of Afc follows from Remark II. II □ 

The regularity of f/j, at A = A^ is assumed throughout this paper. If we 
allow that some Afc coincides with a singularity of Vy^ then the x-dependence 
of Afc in principle can be different from (13.41) and we get an additional freedom. 

3.2 Neugebauer's approach 

A simple but quite general method to construct polynomial Darboux-Backlund 
transformations has been proposed by Neugebauer and his collaborators 
|17l [511 [52] , see also [3ll[60]. We are going to find conditions on polynomial D 
implying that divisors of poles of IJ^, and Vy coincide (compare Remark ll.2p . 
From ^ ,y = Uy'^ we get 

det^(A) detD 

where by D'^ we denote the matrix of cofactors of D. Obviously D'^ is also a 
polynomial in A. 

If Uy are rational functions of A, then given by (13. 7p are rational as well 
(because D and D^^ are rational). Therefore the only candidates for poles 
of Uy are poles of Uy and zeros of detD (i.e., A^). The necessary condition 
for the regularity of Uy ai X = Xk is 

^,y{Xk)^%Xk)=0 . (3.8) 
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If Afc is a simple zero of det-D(A), then the condition fl3.8p is also sufficient. 

Following [17j, we will find another, more constructive, characterization 
of the condition (13.81) . If detD{Xk) = 0, then we have also 

det§(Afc) = (3.9) 

(because ^I^(A) = £)(A)\1/(A)). We assume that the function \1/(A) (known as 
a "background solution" or a "seed solution") is non-degenerate at A = A^. 

As a consequence of (13. 9p . the equation \^(Afc)pfc = has a non-zero solu- 
tion pk e C" (where, in principle, pk can depend on x) . Then, we compute: 

(Afc)Pfc = U^iXk)^iXk)Pk = , 

where we took into account that ^I^(A) satisfies (11.21) . Thus we have: 

^(Afc)pfc = ^„(Afc)pfc = , (3.10) 

what implies (13. 8p . as one can see from the following fact of linear algebra 
(HZ], see also [31]). 

Lemma 3.5. Let us consider two degenerated matrices X and Y . Suppose 
that there exists a vector p such that Xp = and Yp = 0. Then: YX'^ = 0. 

Proof: Let us perform computations in a basis (ci, . . . , e„) such that Ci = p. 
Then all elements of the first column of matrices X, Y are equal to zero. 
Thus, using the definition of the cofactor, we easily see that the rows of Y'^ 
(except the first row) have all entries equal to zero. Hence, XY'^ obviously 
yields zero. □ 

Lemma 3.6. The vector pk such that '${Xk)pk = is defined up to a scalar 
factor. If Xk is a simple zero, then we can choose this multiplier in such a 
way that pk = const . 

Proof: We differentiate the equation defining pk. {Xk)pk + ^{Xk)pk,u = 0. 
Hence, ^{Xk)pk,u = 0, which means that pk,u is proportional to pk (provided 
that Xk is a simple zero of detZ)(A)). Thus pk,u = fkuPk, where fku are some 
scalar functions. From the identity Pk,ufi = Vki^iv it follows that fku^^i = fk^lw 
Therefore, there exists ipk such that fkv = Vkiv Hence, Pke~'^'' does not 
depend on x. □ 
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Corollary 3.7. Polynomial Darhoux matrix liS. 1\) can be constructed as fol- 
lows. In the isospectral case we choose Nn pairwise different complex num- 
bers Xi, X2, . . . , Xnu and Nn constant C"'-vectors pi,P2, ■ ■ ■ ,PNn- We also 
choose the matrix Tq ("normalization matrix"), detTo 7^ 0. Matrix coeffi- 
cients Ti, . . . , T/v are computed from 



where \I/(A) is given ("seed solution"). In the nonisospectral case we choose 



For a fixed k the equation (13. lip consists of n scalar equations. Thus 
we have a system of n'^N equations for N unknown matrices n x n. In the 
generic case such system should have a unique solution. 

The freedom in choosing To corresponds to a gauge transformation. Note 
that an identical situation has place in the case of the binary Darboux matrix, 
where J\f is, in general, undetermined. Usually it is sufficient to put Tq = I 
("canonical normalization"). If this choice leads to a contradiction (i.e., 
the Darboux matrix with the canonical normalization does not exist), then 
we may relax this assumption and search for Darboux matrices with more 
general normalization. 

The case det Tq = can be treated in a similar way but with one excep- 
tion: the total number of zeros is smaller than Nn. As an example of such 
situation we will present elementary Darboux matrices, see Section |H 

3.3 Explicit multisoliton formulas 

Let us introduce the notation 



where (pk G are column vectors. We assume det Tq ^ and denote 



D{Xk)^{Xk)pk = , {k = l,...,nN) , 



(3.11) 




ipk := ^'(Afc)pfc , 



(3.12) 



{j = l,...,N) , 



(3.13) 



where Tj are defined by (13.11) . The equations (13.111) read: 





(3.14) 
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where M = nN. After the transposition we get: 

N 



N-j 
k 



(3.15) 



It is convenient to solve these equations in the matrix form: 



\ Am Vm • • • >^M^M ^ii J 



(3.16) 



Usually, in practical applications, one uses Cramer's rule to express 6'^ in 
terms of determinants, compare [521 EHl CO] ■ 

Having coefficients we can apply the Darboux transformation to Lax 
pairs of prescribed form. As an illustrative example we present the simplest 
but very important case (linear in A): 



Vx = UqX + Ui , 
The equation (II. 4p for v = 1, i.e., JJiD = DUi + D,i, yields: 



(3.17) 



N 



N 



N 



k,l 



(3.18) 



{uoX + u,) A^-'Tfc = A^"'T, (woA + wi) + ^ A^-'^T, 

fc=0 fc=0 fc=0 

Considering coefficients by A^^^ and A^, we get explicit formulas for the 
transformed fields Ui and Uq. 



Uq = TqUqTq \ ui = TqUiTq ^ + [TiTg \ Mo] + ^ci ^ • 



(3.19) 



In the classical AKNS case uq = ia^ = diag(2, —i) and it is sufficient to take 
the canonical normalization Tq = /. Therefore, we get 



Uo = Uo , ui = ui + [Ti, ia-s] , 
where Ti = 9i can be explicitly computed from fl3.16p . compare 



(3.20) 
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4 Elementary Darboux matrix 

The elementary Darboux matrix is linear in A and its determinant has just a 
single simple zero. This case is mentioned by Its [35] and discussed in more 
detail in, for instance, pU |38]. An obvious way to produce matrices of this 
type is to take matrices with a single entry linear in A and all other entries 
A-independent. In this paper we confine ourselves to elementary Darboux 
matrices for n = 2. They can be represented in the form 



where AT, M. do not depend on A. As a simple exercise (compare Corol- 
lary 13. 7p we can express the coefficient a by \I' evaluated at Ai, namely: 



where pi is a constant vector. 

4.1 Binary Darboux matrix as a superposition of ele- 
mentary transformations 

Theorem 4.1. In the case n = 2 any binary Darboux transformation is a 
superposition of two elementary Darboux transformations. 

Proof: We will show that 



is a binary Darboux matrix (A/i,A/2,A/l are non-degenerate matrices which 
do not depend on A). First, performing the multiplication in (14.31) . we get 




(4.1) 




(4.2) 




(4.3) 



D=Ar(A-Ai + (Ai-A2)P) 



(4.4) 



where 




(4.5) 
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and AA = Ai — A2. Then, we easily check that = P. 

The coefficients a,P can be expressed by ^ evaluated at Ai, A2. Indeed, 
denoting 

^(A.K =( i'] , (4.6) 



and using equations (13. lip , we obtain 

Vi n 66AA 
a = ^ , f^= 7 7- ■ (4-7) 

The projector P reads 

P = ^-—M-' ( "^'^^ \ M . (4.8) 

If = J, then ^(Ai)pi G kerP and '^{\2)P2 ^ imP. Therefore, the bi- 
nary Darboux matrix with P given by (14. 8 p is a superposition of elementary 
transformations (14 .Sp with M. = I and a, (3 given by (14.70 . □ 



4.2 KdV equation 

The Darboux transformation for the famous Korteweg-de Vries equation is 
almost always presented in the scalar case, see [16]. The matrix approach 
is less convenient. However, having in mind a pedagogical motivation, we 
are going to show in detail that the matrix construction works also in that 
case. It is interesting, that in this paper we do not need the "KdV reality 
condition" (usually used in earlier papers, compare [HI [271 176]). 

The standard scalar Lax pair for KdV equation consists of the Sturm- 
Liouville-Schrodinger spectral problem and the second equation defining the 
time evolution of the wave function: 

-ip^ii+uip = Xip , ?/',2 = -4?/',iii +6u?/',i +3m,i . (4.9) 

The compatibility conditions ip,ii2 = V',211 yield the KdV equation 

u,2 -6mm, 1 = . (4.10) 
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The Lax pair fl4.9p can be transformed, in a standard way, to the matrix 
form 



'2 1 -^AX^ + 2u\ + 2u^ - u,n u, ' 



(4.11) 



where 



^'=(v'J' ^'=(1)' 

and V^, are hnearly independent solutions of (14.91) . 
Lemma 4.2. Suppose that 

where u,a,b,c do not depend on X. Then, the compatibility conditions 
U,2 —V,i + [U, V] = uniquely yield: 

a = u,i , b = 2u , c = 2v^ — . (4-14) 



i.e., U, V given by I!i4-1^ ) are identical with the Lax pair ( [^.ii| ) for the KdV 
equation. 

Proof is straightforward: compatibihty conditions reduce to (14.101) and (I4.14p . 

4.3 Elementary Darboux matrix and the classical Dar- 
boux transformation 

We will compute the action of the elementary Darboux transformation in the 
KdV case, compare [27] • We assume 

D^m{'-_^^\) (4.1B) 

where N (detA/" ^ 0) does not depend on A and a is a function to be 
expressed by ^'(Ai), namely 

D(Ai)^(Ai)pi = , (4.16) 
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where pi is a constant vector. We denote 

= ^(Ai)pi (4.17) 
The constraint (14.161) (with D given by (14.1 5p ) is equivalent to 



a = — = ^, (4.18) 

where t/'i satisfies (14. 9 p with A = Ai (i.e., ^/'i is a hnear combination of ipi 
and 01 ). The function a satisfies the following system of Riccati equations: 

a,i = M — Ai — , 

(4.19) 

a,2 = (2^2 - +2nAi - 4Ai) + 2u,i a - {2u + 4Ai)a^ , 

which can be obtained directly from (14. lip . 

The elementary Darboux transformation for U, V (i.e., the formulas (11.40 
with D given by (I4.15P ) reads: 

^-x^-^(o iy'-^{7 1')^-'^^"-^"- 

(4.20) 

where 

Ml = (u - Ai - - a,i ) AT J |j ^ A/"-! , 

Ms = (-a,2 +c + 2aa - ba^ - iXj + 2Ai(m - 2a^)) A/" J [J ^ TV"^ , 

(4.21) 

and Vq does not depend on A (its explicit form follows from Lemma [4.21 and, 
therefore, is automatically preserved by the Darboux transformation). 

The necessary condition for the Darboux transformation is vanishing of 
residua Mi,M2 (what is equivalent to (I4.16P and, as a consequence, to the 
Riccati equations (14.190 ). 
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In order to assure the Darboux invariance of the coefficients by A in [/ 
and by in V we have to impose some constraints on the normahzation 
matrix A/" (compare |41j), namely 

^(r.)^--( (4.22) 



what imphes the following form of Af: 

^ = / ( _°1 J ' (4-23) 
where /, 7 are functions of x. Now, the transformation fl4.2UI) becomes 



\u-X a-j J ^ f \0 1 

where 

u = \i + 27a — 7^ — ^^-^ 
V = Xi + 27a; — 7^ — 26 



(4.24) 



(4.25) 



Comparing (14.241) with (14.131) we find the remaining constraints on the form 
of the Darboux matrix: 

/ = const , 7 = 0;, 2v = u . (4.26) 

We assume / = — 1. By virtue of (14.141) b = 2u, and we easily verify that the 
constraint 2v = u coincides with the first Riccati equation (14.191) . 

Corollary 4.3. The elementary Darboux matrix for the KdV equation is 
given by 

where a is computed from (^T7^, see also (^T7^. 
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The transformation of u can be obtained from fl4.19l) and fl4.25p . see 
fl4.28p . Taking into account fl4.12l) we get the transformation for ip. 

Corollary 4.4. The elementary Darhoux matrix \4-21^ generates the classi- 
cal Darhoux transformation: 

(4.28) 

M = M — 2a,i = -u — 2(ln-?/;i),ii , 
where -ipi = ip{x, Ai) satisfies ^4-^ - 

Formulas fl4.28p were first obtained by Gaston Darboux [23], see also |46j . 

Proposition 4.5. D = -Dq,,Ai given by \4-21^ has the following properties: 
D~\^ is equivalent to D_a,Xi and Dp^\^Da,\^ = M{\ — Ai + M), where 
M '= (Ai - A2)P (and = P) for A2 ^ \i and = for A2 = Ai. 



Proof: by straightforward computation. First, D-q^Ai = Then, 

-1 

a + f3 



DpMDaM = { °i ) (A - Ai + M) (4.29) 



where 

M-( -«(« + /^) « + A .43(.^ 

~ V «(Ai - A2) - a\a + P) «(« + /?) - (Ai - A2) J' ^ ' 

and we easily verify that = (Ai — A2)M, which means that (for A2 7^ Ai) 
M = (Ai - A2)P (where P^ = P), compare □ 



4.4 Nilpotent Darboux matrix and classical binary Dar- 
boux transformation 

Let us consider the Darboux matrix of the form (12.31) . In the case n = 2 the 
nilpotent matrix M (M^ = 0) can be parameterized as 

M = ( I ) (4.31) 

where g, a are some functions. 
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Considering the transformation (ll.4p we have to demand that 11^, are 
regular at Ai by canceUing the pole of second order at A = Ai. We get two 
conditions: 



M,, +[M, f/,(Ai)] - MUl{\i)M = , 

(4.32) 

M,^M + Mf/^(Ai)M = , 

where the prime denotes differentiation with respect to A. The second set of 
equations turns out to be a consequence of the first equations (it is enough 
to multiply them by M from the right). 

The classical binary Darboux transformation is usually defined only for 
the time-independent spectral problem [IB]. Therefore, in order to show 
that the considered transformation (12.31) coincides with the classical binary 
transformation it is sufficient to confine ourselves to u = 1. The equations 
(I4.32P (for u = 1) can be rewriten in terms of g, a: 

g,i -2(T^ + ^^ = , 

g,i a + ga,i +ag^ - g{u - Ai + a^) = , (4.33) 

g,i + 2gaa,i +a'^g'^ - 2ga{u - Ai) = , 
Using the first equation we can reduce the last two equations to: 

a,i+(7^ -M + Ai = . (4.34) 
Therefore 



a = 




(4.35) 



where ipi satisfies the first equation of (14.91) for A = Ai, compare (I4.18p . 
(I4.19p . Taking into account (I4.32p we rewrite (ll.4p for f/i = f/ = uqX + ui as 

uq = MuqM-^ , ui=J\f,iM-'^ +M{ui + [M,uo])M-'^ . (4.36) 

In the KdV case, see (I4.13p . the first equation of (I4.36P is satisfied for 




(4.37) 
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Then, the second equation of fl4.36p reduces to: 

■u = M + 7,1 +2(77 + 7^ , g = -'y . (4.38) 
Taking into account the first equation of fl4.33p we finally get get 

u = M + 27,i , 7,1 -2(T7 - 7^ = , (4.39) 
where a is given by fl4.35p . Therefore, the last equation is equivalent to: 



9a; Y 7 ^ 
which means that 



(4.40) 



7 



Co- I^f 



u = u 



dx' 



In 



Co 



(4.41) 



where cq is a constant of integration. The last formula coincides with the 
classical binary Darboux transformation for the Sturm-Liouville-Schrodinger 
spectral problem 



Corollary 4.6. The nilpotent Darboux matrix Ii2.3\) generates the classical 
binary Darboux transformation. 

The "second" binary Darboux transformation, introduced in [83j, corre- 
sponds to the choice cq = 1. 



5 Fractional form of the Darboux matrix 

Another popular representation of the Darboux matrix (with nondegenerate 
normalization) is decomposition into partial fractions [56l [79| [80l [82] : 

D=Af{I+-^ + ...+ ^ 



D-' = {I+-^ + ... + -^)Af-' . 

A — /ii A — /i^r 

In principle the numbers of poles of D and could be different but here, 
following other papers, we assume the "symmetric" case fl5.ip . 
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We will denote by the Darboux matrix in the fractional form with the 
canonical normalization (in other words, D = MDq). The form fl5.ll) of D 
and imposes restrictions on Ak and Bk implied by equations DD^^ = I 
and D-^D = /, see 

Multiplying D by the least common multiple of the denominators we 
obtain the equivalent polynomial form D{X) (a polynomial of A^th degree). 
The determinant detD{X) is a polynomial of degree Nn vanishing at poles 
of D and D~^, i.e., at A = and X = fik {k = 1,...,N). The sum of 
multiphcities of all zeros of detD{X) equals Nn. Therefore, for rz = 2 all 
zeros are simple, while for n > 2 some of them have to be multiple zeros. 

The fractional form is convenient in the case of some reductions (e.g., 
orthogonal or unitary), when the eigenvalues A^ {k = 1, . . . ,nN) can be 
naturally divided into pairs Xk, fJ^k- 

5.1 Zakharov-Mikhailov's approach 

We start from fractional representation of the Darboux matrix (15.11) . where 
Ak, Bk have to satisfy constraints resulting from the condition DD~^ = /: 




{k = 1,...,N). We assume the nonisospectral case and demand that Ui, 
defined by (ll.4p have the same form as U^. In particular, it means that 
the right-hand sides of (ll.4p have no poles. Equating to zero the residua at 
A = Xj and at A = fik, we get 

(A,,^ +A,t/(A,)) [i + Y: + {L.{\) - A,,. ) E ' = 0, 

^ A \ ^ A B 
I + Y] V MBk - Bk,,, ) - (L^(^fc) - fik,u ) Y] 7 ^^T^ = 0, 

(5.3) 
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for j,k = 1, . . . , N . Multiplying first equations by Aj from the right and the 
second equations by from the left, and then using (15.21) . we obtain 



(5.4) 

which is satisfied when (13.41) (and similar equations for fik) hold. Note that 
we derived here a proposition analogical to Theorem 13.41 

In order to solve the system (15. 2p . (15. 3p we assume (13.41) and represent 
Ak, Bk as follows: 

Ak = \sk){ak\, Bk = \bk){qk\ (5.5) 

where \sk), \bk ) are matrices built of linearly independent n-component 
column vectors and { qk\, { CLk\ are matrices built of linearly independent 
n-component row vectors. In other words, all these matrices have maximal 
rank. In particular, | Sk ) and ( ak \ have the same rank (denoted by ikAk) 
but (in general) different than the rank of 1 6^ ) and ( qk \ (denoted by ikBk). 
Using the notation (15. 5p we rewrite equations (15. 2p and (15. 3p as follows: 



\sk){ ak\DQ^{\k) = , Do{nk)\bk ){Pk\=0 , 
DQ^{Xk)\sk )( afcl = , \bk ){Pk\Do{lJ'k) = , 

Sk ),u {ak\ + \sk )( ak\,u+\sk )( ak\Uu{Xk)^DQ^{Xk) = 



Do{fik)(^- \bk ),u {qk\ - \bk )( qk\,u+\bk )( qk\Uu{fik)^ 
where k = 1 N and u = 1 m. Moreover, 



(5.6) 



(5.7) 







Do{,k) = I / + E , ^o-^(A.) =(l + t ) • (5-8) 

* ~1 ^'k-XJ J \ ~t ^k-fij ' 



Lemma 5.1. If \ a ) and { b\ have the maximal rank, then: 

|a)(6|=0 \a) = or (6|=0. (5.9) 
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Proof: immediately follows from the definition of the maximal rank. All 
columns of |a ) (and all rows of ( b\) have to be linearly indpendent. □ 

Using (15.61) and applying Lemma [STT] to equations (15.71) . we get the fol- 
lowing linear system: 

{ ak\,u= -{ ak\Uu{Xk) , \hk),u=Uy{iXk)\hk) , (5.10) 
which is satisfied by: 



{ako\^~\\k), l&fe ) = *(/Ufe)|6feo ) , (5.11) 



where ( a^o \ and | h^Q ) are constant. If \E' is regular at and /i^, then the 
solution given by (15.111) is general (compare Section [23]) . 

5.2 Symmetric representation of the Darboux matrix 

We proceed to derive compact formulas for the remaining ingredients of 
namely for ( \ and \pk)- Taking into account Lemma 15.11 we can simplify 
equations (15. 6p : 

N N 

(cik' ' 



+ ^Mfc,-( g,| =0 , \bk)-J2\^j)^Jk = 
j=i j=i 

N N 

{ qk \ + '^Kkj{ ak \ = , \sk) - ^\bj )Kjk = . 



(5.12) 



where 

A4, = i5iiM. A> = Mii). (5.13) 

Afc ~ fij fij — Ak 

The expression (a^ | bj) denotes matrix multiplication: (a^ | bj) = { ak\\bj ) 
(for any fixed j,k). The resulting matrix is not necessarily quadratic. The 
number of its columns is Tk{Bj) and the number of its rows is rk(y4fc). Simi- 
larly, {qj I Sk) is also a matrix (for any fixed j, k). 

Remark 5.2. Matrices Mjk form the so called "soliton correlation matrix" 
M which has XljLi ^k(-Bj) columns and Yl!k=i'^^{^k) rows. From \5.1^) it 
follows that 

k = M-^ . (5.14) 
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Therefore, M and K have to he quadratic matrices, i. e. 



N N 



^rk(A,) = J]rk(S,) . 



(5.15) 



k=l 3=1 



The soliton correlation matrix M is a Cauchy-like matrix (compare [53] ) 
which has been reobtained several times in various particular cases (see, for 
instance, [121 [321 [621 [76] ) . 

Corollary 5.3. The symmetric form of the multipole Darhoux matrix is given 
by 



where M is a normalization matrix (we assume detTV T^Oj, K = M ^ , M 
is given by ^5.13\) . and \ hj ) , { aj \ (j = 1, . . . , N) are expressed by Ii5.11\) . 

5.3 How to represent A^-soliton surfaces? 

Iterated Darboux matrix is a composition of binary Darboux transforma- 
tions (see, for instance, [lOl [56] ): 




(5.16) 




D=J\f{l + 




Pn 



) 




(5.17) 



where projectors Pk are defined by 



kerPfc = ^fc-i(Afc) , imP^ = ^fe„i(/ifc) , 



(5.18) 



where are defined by: \E'o(x, A) = \E'(x, A) and (for k ^ 1): 




(5.19) 
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In this case f ll.6p yields 



N 



^ (A - Afc)(A - /ifc) 



(5.20) 



Note that the formula (15.201) does not contain Af. Indeed, gauge equivalent 
linear problems have identical soliton surfaces, see [73] . 

Remark 5.4. The determinant of the iterated Darboux matrix ( [5. i7| ) can be 
easily computed (compare fl4^): 



N 



detD = detATjJ 



A -A. 



dim imPfe 



(5.21) 



Usually the formula (15.201) is used in the isospectral SU{n) case, when 
A,^ = 1, fik = Xk (see Section [6]), and TrF = TrF = (what can be attained 
by multiplying D by an appropriate factor /, see (11.71) ): 



lA-A 



L.-^^Wi('^}^I-P,)^,_,W, (5.22) 

k=l r- --kl \ n / 

(in this case projectors are orthogonal, = Pk, see for example [HI [391 [71]). 
The sum on the right-hand side of (I5.22p consists of traceless components of 
constant length (using the Killing-Cartan form a ■ b = — nTr(a6) as a scalar 
product in su{n)), see [H]. Thus the formula (I5.22p generalizes the classical 
Bianchi-Lie transformation for pseudospherical surfaces [72], [71]. 

The formula (I5.20p is not manifestly symmetric with respect to permuta- 
tions of Afc. The symmetric formula for the Darboux-Backlund transforma- 
tion for soliton surfaces can be obtained by substituting (15.161) into (11.60 . 

Theorem 5.5. The symmetric representation for N -soliton surfaces has the 
form: 



(the notation is explained in Corollaru \5.3\ see also Section 1.4\ ) 
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Proof. We compute D ^D,x where D is given by fl5.16p : 

J,k=l ^ ' t,j,k,l=l 

We use f l5.13p and perform the summation over i,l in the second component: 

N N N 

J2 Kji{\,-fii)M,iKik = Y,K,X5tk-Y,SfiiiiKik = K,k{\k-fij), (5.24) 

i,l=l i=l 1=1 

where 6{l, 5^ are natural generahzations of Kronecker's deha (e.g., 6^j^ is unit 
matrix of order ik{Ak) and Sfj is unit matrix of order Tk{Bj)). Finally, by 
virtue of an obvious identity (A — fij) — (A^ — fij) = A — Afc, we get (I5.23p . □ 

The expression (15.230 is a generalization of the symmetric formulas for 
A^-soliton surfaces which has been earlier obtained in the S'u(2)-AKNS case 
see also 



Proposition 5.6. The Darhoux matrix fD, with D given by ( 15. and f 
given by 



f 



friX-XkY^ 



transforms traceless F into traceless F. What is more, det{fD) = det A^. 
Proof: If D is given by (15.161) . then, using (I5.23p . we compute: 

A,c Vj>^i(A-Afe)(A-/i,)J A. A -Afc A-/ifc 

(5.26) 

where we took into account that X^jLi ^kjKjk is the unit matrix of order 
rk(ylfc) and ^^Li -^ifc^fcj is the unit matrix of order ik{Bj). Multiplying D 
by a A-dependent function we can change F — F, by virtue of (II. 7p . In order 
to get TrF = TrF, we have to take / such that the right-hand side of ( 15.26^ 
equals n(ln/),A. Hence we get (I5.25p . 
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Surprisingly enough, in this way we can compute also the determinant 
of D given by flETBD . Indeed, from we have TT{{fD)-^{fD),i;) = 
(provided that TtF = TtF) and then Theorem 13.21 implies that det{fD) 
does not depend on ( (and is A-independent, as well). Therefore we can 
evaluate det(/Z)) at A = oo. Thus we obtain det{fD) = detM. □ 

Let \1' = fD"^, where D is given by fl5.16p and / is given by (15.251) . For 
simplicity we assume also rk(Afc) = rk(i?fc) = r^. Then: 

F = F + X,,^~\X) y ( - - 1^^- )KM a,\ \ ^^^^ ^^27) 

(A-Afc)(A-^j) J 

where 6jk is Kronecker's delta. Note that the symmetric form of (I5.22p is 
given by the specialization of the formula (I5.27P to the case /i^ = Xk- 

Another representation for multisoliton surfaces can be derived from the 
polynomial representation of D. In order to compute A^-soliton addition to 
the surface F := \I'~^\I',;^ |a=Ao assume the Darboux matrix in a general 
form 

N 

D = ^Tfc(A-Ao)'= . 

k=0 

The matrices Ti, . . . , T^r are computed from the following linear system: 

N 

Tk{X, - Ao)'^(A,K = , (z/ = 1, . . . , iVn) , 

A;=0 

where X^ E C and Pi, G C" are constant, and Tq is a given normalization 
matrix. Of course, one should take care of reductions what can result in 
some constraints on Aj,, Pi, and also on Tq, see Section [61 The formula (ll.6p 
assumes the form: 

F = F + ^-\X) T^^Ti ^(A) = F + ^-^(A) Oi ^(A) , (5.28) 

where 9i is given by (I3.16p . Note that also in this case F does not depend on 
the normalization matrix To (a change of Tq implies such change of Ti that 
6*1 remains unchanged, compare Section [3l3|) . 
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6 Group reductions 



The so called reduction group was introduced by Mikhailov [49j and detailed 
description of various reductions is given in [50], [80], see also [l^. Group 
reductions (under a different name) found a rigorous treatment in the frame- 
work of the loop group theory [291 EHl EE] , compare Section I7.3[ 

In this section we describe several important types of reduction groups. 
We consider only the case of non-degenerate normalization det A/" 7^ (which 
means that detD{X) is a polynomial of degree Nn). The most convenient 
form of the Darboux matrix depends on the reduction. The polynomial form 
f l3.1l) is very good for reductions to twisted groups, the fractional form (15. ip 
(and especially its symmetric version (15.161) ) is appropriate for unitary and 
orthogonal reductions. Then, as an example, we present in more detail the 
principal chiral model (sigma model) and its reductions. The symmetric form 
(I5.16P is of great advantage in this case. 

6.1 Reductions to twisted loop groups 

Twisted loop groups are defined by \l'(ci;A) = Q'^{X)Q^^ where ou = exp ^ 
(hence = 1) and, necessarily, = I (we assume also Q = const), 
compare [21]. An important example, two dimensional Toda chain (then 
i^' = n), is discussed in detail in [50], using the fractional representation 
of D (15. ip . Here we present a different approach, based on the polynomial 
representation. 

Usually it is better to consider some natural extensions of loop groups 
which follow from the form of the linear problem. The starting point is the 
assumption about the form of the linear problem (i.e., U^, are constrained to 
the corresponding Lie algebra): 



which implies (^(u;A))„,= Uu{ljX)'^{ujX) = QUu{X)Q''^'^{ujX). Hence 
{Q-'^iuX)) U,{X) {Q-'^iouX)) , 



which means (see Remark II. 4p that Q ^\l'(a;A) = \l/(A)Co(A), where the 
matrix Co (A) does not depend on x. Therefore 



U,{u;X) = QU,{X)Q 



(6.1) 



^(^A)=g^(A)Co(A) , 



(6.2) 
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and similar equation for \^ = Z)\l/ (with a different Co{X), in general). There- 
fore: D{ujX)Q<il{X)CoiX) = QD{X)^iX)Co{X). In order to eliminate ^(A) we 
have to assume that Co(A) = 7o(A)Co(A), where 70 : A ^ 7o(A) G C is a 
rational complex function of A. Then 

b{ioX) = 7o(A)Q^(A)g-i ^ 

(6-3) 

detD{uX) = (70(A))" detD(A) . 

Remark 6.1. Computing D^u'^X) ,. . . ,D{uj^ X) we obtain a necessary con- 
straint for 7o .■ 

7o(A)7o(c^A)...7o(^'^"'A) = 1 . (6.4) 

This constraint is satisfied by any meromorphic function such t/iat 70(00) = 1 
and all its zeros and poles coincide with some zeros of det D{X) . Note that 
the matrix Co(A) also is not arbitrary but satisfies an analogical constraint. 

We make usual assumptions: 70(A) = 1 and Co(A) = Q^^ (then \I' and 
D are fixed points of the reduction group [50j, or, in other words, D and \& 
take values in the loop group). Then 

^{u^X) = g^*(A)Q-*= , D{u;'X) = Q^b{X)Q~^ , (6.5) 

for A; = 1. 

Lemma 6.2. Let 70(A) = 1. // detZ)(Ai) = 0, then det £)(cu'=Ai) = for 
k = 1, . . . , K . Multiplicities of all these K zeros are identical. 

Therefore, if -D(A) is a polynomial of order (and, as a consequence, 
detD{X) has the order Nn), then Nn has to be divided by K, i.e., there 
exists an integer N such that Nn = NK (in the two-dimensional Toda chain 
case N = N). Moreover, fl6.5p (evaluated at A; = 1) imply that 

cu^ToQ = QTo , (6.6) 

where To is the normalization matrix, compare (13. ip . We have to demand 
that this equation has a solution To 7^ (otherwise, the Darboux matrix 
cannot be a polynomial of order N). Certainly (16.61) has a solution for 
such that uj^ = 1 (in fact, this assumption was done in |50j). 
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Corollary 6.3. //70(A) = 1, then the set of zeros of det-D(A) is given by 
{iu'^Xjl = 0, 1, . . . ,K - 1; j = 1, . . . ,iV}, where Xj G C. 

Equations (13.111) . defining the Darboux matrix, can be rewritten as fol- 
lows (taking into account (16.51) ): 

= D{uj'Xj)^i/{uj''Xj)pjk = Q''D{X,)<i/{X,)Q-%k (6.7) 

For simplicity we assume the generic case, i.e., all zeros cu'^Aj are pairwise 
different (and, as a consequence, simple). Then the kernels of D{Xk) are 
one-dimensional, which means that '^{Xj)Q~^pjk is proportional to '^{Xj)pjo. 

Theorem 6.4. Assuming that \6.0^) has a solution Tq ^ we construct 
the Darboux matrix (a X-polynomial of order N) according to Corollary \3. 7| 
taking into account that its zeros are given by oj^Xj (see Corollary \6.3\) and 
the corresponding eigenvectors are related bypjk = Q'^Pjo (k = 0,1, . . . , K—1; 
j = 1,...,N, where N = Nn/K). This Darboux matrix preserves twisted 
loop group constraints lid. 1\) and Ii6. 5\) . i.e., UyiuX) = QU^{X)Q~^ , etc. 

Remark 6.5. In the noniso spectral case twisted reductions impose constraints 
on the form of Ly. If \3.Ji\) are satisfied, then also ujXk,u = Lu{x,u}Xk) ■ Hence 
we get the constraint: ujLy{x,X) = Ly{x,ujX). 

The particular case K = 2 (i.e., uj = —1) is very popular (e.g., this 
reduction is necessary to derive the standard linear problem for the famous 
sine- Gordon equation [56l[60], see also \18\)- This case can be generalized by 
admitting a A-dependence of Q (actually such generalization can be done for 
any K but the results have more complicated form, so we omit them). One 
can easily see that Q = Q{X) has to satisfy 

g(-A)g(A)=^o(A)/, (6.8) 

where {}o is a scalar function such that 'i9o(— A) = 'do{X) (in particular, we 
can take ^o{X) = 1). Assuming 70 = 1 we have detD{—X) = detD(A) 
which means that zeros of detD{X) appear in pairs A^' = —A/,.. Constant 
eigenvectors pk' and pk satisfy ^{Xk)pk = and '^{Xk')pk' = which implies 
'^{Xk)Q{—Xk)pk' = 0. If the zero A^ is simple then pk' = Q{Xk)pk (the 
eigenvectors are defined up to a scalar constant factor, therefore we omitted 
the factor ^o{Xk)). Moreover, the condition (16. 6p should be replaced by 
uj^TqQoo = QooTq, where Qoo is either Q{oo) or the coefficient by the highest 
power of A in the asymptotic expansion of Q{X) for A —>■ 00. 
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6.2 Reality condition 



The condition 



f/,(A) = U,{\) (6.9) 

(where the bars denote complex conjugates) simply means that all coefficients 
of matrices Uy are real. Considering (11. 2p we get 



^(A) = ^(A)C(A) (6.10) 



where C(A)C(A) = /. Applyi.ng fl6A0l) to ^(A) = D{\)^!{\) we obtain the 
corresponding constraint on D 



DCX)=i{X)D{X) , (6.11) 



where 7(A) is a scalar rational function which satisfies 7(A)7(A) = 1. Hence 



which means that 7 = w{X)/w{X), where ti'(A) is a polynomial of degree 
K ^ N (provided that D is a polynomial of degree N). Then detD{X) has 
K zeros of multiplicity n and {N — K)n simple zeros. The set of simple zeros 
is invariant with respect to the complex conjugation, i.e., they are either real 
or form pairs of conjugate numbers. 

We assume the simplest case 7(A) = 1 (and also C(A) = C(A) = 1). Then 
all zeros of detD{X) are simple, and either Xj G M (then pj = pj) or there 
are pairs A^/ = A^ (then pk' = Pk)- 



6.3 Unitary reductions 

Unitary reductions (which sometimes are also referred to as reality condi- 
tions, see for instance [70]) are defined by 

UlCX) = -HU,{X)H~^ , (6.13) 

where the dagger denotes the Hermitean conjugate, if is a constant Her- 
mitean matrix (if^ = H) and v means complex conjugate (necessary if x^^x^ 
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are complex, e.g., usually = z, x"^ = z in the case of cliiral models, 
discussed in Section [n3D- Using fl6.13l) we obtain from (11.21) : 

(*t(A)),,^ (vl/(A)„)t = vl/t(A)f/t(A) ^ -^^i\)Hi\)U^{\)H~\\). (6.14) 

Taking into account the well known formula for differentiating the inverse 
matrix (i.e., (\l/~^),i,= — \l/~^\E',j, we transform fl6.14p into 

H-\{^\\))-'),, = U.{\)H-\¥C\))-' , (6.15) 

and, comparing fl6.15p with (11.21) . we get 

(vI/t(A))-i = i7vi/(A)Co(A) , (6.16) 

where ^I/(A) solves the system (II. 2p and Co(A) is an x-independent matrix. 
From (I6.16P we can derive Cq(A) = Co (A). ^(A) satisfies the constraint (16.161) 
with Co in the place of Co. Assuming Co(A) = fco(A)Co(A), where ko{X) is an 
x-independent scalar function, we derive the condition 

Z)t(A) = ko{X)Hb'\X)H-^ , (6.17) 

which is necessary for D to be the Darboux matrix in the case of unitary 
reductions. We point out that A;o(A) has to be a rational function. The 
simplest choice /co(A) = 1 is not possible. Indeed, from (I6.17P we obtain 

(fco(A))" = detD{X) detD{X) . (6.18) 

Hence /co(A) is a polynomial of degree 2N (provided that detTo 7^ 0). More- 
over, fco(A) = fco(A) which means that fco(A) is a polynomial with real coef- 
ficients. The set of its zeros is symmetric with respect to the real axis. We 
confine ourselves to the case of ko{X) without real roots, i.e. 

A;o(A) = I det To|2(A - Ai)(A - Ai) ... (A - Ajv)(A - A,v) 

where is the degree of the polynom -D(A). Then 

detD = (detTo) (A - Ai)"-''^^ - ^i)''^ ... (A - A;v)""''^(A - Ajv)""^ , 

where dk are some integers, 1 ^ ^ n — 1. Note that for n > 2 the zeros of 
detD rule, degenerated. If A^ are pairwise different, then dividing 

D by (A — Ai) . . . (A — Aat) we obtain the matrix D (equivalent to D) bounded 
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for A — s> oo (liiiiA^oo D{\) = Tq), with singularities at A = (A; = 1, . . . , A^). 
Taking into account detTo 7^ 0, we get: 



where Af = Tq and are some matrices dependent on x. The inverse matrix 
D^^ has poles at A = A^ {k = 1, . . . , N). Therefore D is exactly of the form 
(15.11) with /ifc = Afc. Note that (16.171) can be rewritten as: 



In what follows we use natural notation: {a^^ ) := (( a\y, { := {\b ))^ 

Theorem 6.6. The Darhoux matrix of the form ( (5. satisfying the addi- 
tional constraints 



/ifc = Afc, |6,o) = Co(A,)|4o), M^HU = H, rk^ = rk5,, (6.21) 



(k = 1, . . . ,N ), preserves the unitary reduction defined by lid.l!^) and l{6.16\) . 

Proof: We will show that the constraints (16.211) . imposed on D given by 
(I5.16p . are sufficient to satisfy the equation (16.201) . The condition Hk = Xk 
is already assumed. Equating normalization matrices in ( ]6.20p we obtain 
A/"^^ = H^^M'^H . The most convenient way to proceed further is to use 
symmetric form of the Darboux matrix (15.161) . Equating residua at both 
sides of (16.201) we get Bj. = A\ (compare (15.11) ). i.e., 

Af N 




(6.19) 



D-\\)=H-^D\X)H . 



(6.20) 




(6.22) 



In order to satisfy this equation it is sufficient to require 



al) = H\b,) , 



(6.23) 



what implies {b\ \ = { ak\H ^. Indeed, using (15.131) we get Mj^ 



-Ml 



kji 
□ 



and, as a consequence, Kjj^ = —Kkj, compare (I5.14p . 
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Remark 6.7. Assuming Cq{X) = we rewrite the constraint Ii6.16\) as 
\E'^(A)iJ\E'(A) = H, i.e., ^'(A) takes values in the same loop group as -D(A), 
compare Ii6.20\) . This assumption is not very restrictive. It is sufficient to 
impose it on initial data (at x = xq). Then it holds for any x. 

Remark 6.8. In the non-isospectral case the unitary reduction imposes con- 
straints L^{x,X) = L,y{x,X) on the evolution of X, compare Remark \6.5[ 

Remark 6.9. By virtue of Proposition [37g| the Darboux-Bdcklund transfor- 
mation for the reduction SU{n) is generated by the Darboux matrix fD, 
where D is given by Theorem \6.6[ and f is given by liS. 25\) with ikBk = ikAk 
and /ifc = Xk, 

We point out that the case of A;o(A) with real roots is more difficult. 
Usually, the assumption is made that Xk are not real, compare [271 EOl ES] . 
The case of real Ai,/ii is solved and discussed in the case of the binary 
Darboux matrix (A^ = 1), see [Mj. By iterations one can obtain more general 
solutions. However, it would be interesting to obtain a compact form the 
Darboux matrix corresponding to arbitrary set of eigenvalues (symmetric 
with respect to the real axis). 



6.4 Chiral fields or harmonic maps 

As an illustrative example, we will consider the equation 

(<l>,i<l>-i),2+($,2$~'),i = , (6.24) 

which describes harmonic maps on Lie groups (provided that $ assumes 
values in a Lie group G) [76l |77j or, in physical context, principal chiral 
fields [301 US]- Adding another constraint, <l>^ = J, we get chiral fields (or 
sigma models) on symmetric (or Grassmann) spaces [H ED Ell EH] • 

The chiral model fl6.24p is integrable and the associated isospectral Lax 
pair \E',j^ = 11^'^ is of the form 



^,i=Y^^, vl/,2=-^vl>. (6.25) 

The Lax pairs considered in [27] and [77| are equivalent to fl6.25p modulo 
simple transformations of the parameter A. It is convenient to denote 

$(a;) = ^(a;,0) . (6.26) 
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Then A^, = $ ^ or, in other words, 
and the compatibihty conditions, 

^1,2+^2,1=0, Ai,2-A2,l+[Ai,A2] =0 , (6.28) 

rewritten in terms of $ become identical with (16.241) . 

The Darboux-Backlund transformation for $ (in the case of the principal 
GL{n, C) chiral model, where there are no restrictions on $ except non- 
degeneracy) is given by 

$ = L)(0)<l> (6.29) 

where -D(A) is represented, for instance, by the symmetric formula (I5.16p . 

The U{n) reduction is defined by the constraint = J and adding 
det $ = 1 we get SU{n) principal sigma model, see for instance [301 EZ]- 
These constraints are preserved by an appropriately modified Darboux ma- 
trix, see Theorem 16.61 and Remark 16. 9[ 

Chiral models on Grassmann spaces can be characterized by the addi- 
tional constraint: = /. This is a quite non-trivial reduction, worthwhile 
to be considered in detail. 

Proposition 6.10. The Lax pair 1(6. 25\) satisfies the constraints 

U,iX-')=^,,<^-^ + m,iX)^-\ (z/ = l,2) (6.30) 
if and only if $^ = const. 

Proof is straightforward. We check that ^U^^^^ = —Uy iff $^ = const. 
Then, we compute 

^-(^'') = XTFiF ^ "^'^ ~ ^"^^^ ^ "^^^^^^"^'^ ' 

what ends the proof. □ 

The right-hand side of fl6.30p has the form of a gauge transformation. 
Hence, we immediately have the following conclusions. 
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Corollary 6.11. // $2 = I, then 

= f/,(A)$-iM/(A-i) , (6.31) 

which means that 

= <l>^(A)5o(A) , (6.32) 
where S'o(A) is a constant matrix. One can check that So{X~^) = Sq^{X). 
Proposition 6.12. The Darhoux transformation preserves ( fg. 3(A) if 

D{X-^) = <Id(A)$"^ , $ = D(0)$ . (6.33) 

Proof. The constraint fl6.32p for \I' = D"^ reads 

D{X-^)^{X~^) = <Id(A)^(A)So(A) . 
Using we obtain ^M>- Finally, we apply flOB]) . □ 

Corollary 6.13. The formula ^6. 33) implies that divisor of poles of D{X~^) 
has to he exactly the same as divisor of poles of D{X). Inverting lid. 33) we get 
that divisors of poles of D~^{X~^) and D~^{X) also should coincide. Therefore 
both sets of poles, i.e., {Ai,...,AAr} and {fii, . . . , fi^} , o^re invariant with 
respect to the inversion X — > A~^. 

Theorem 6.14. We assume that poles and zeros (XkjfJ^k) of the Darhoux 
matrix ^5.10) can he comhined in the following pairs 

Xk' = , fJ'k' = — , (6.34) 

and A| 7^ 1, /i| 7^ 1. We assume also M = I and 

{ afo \ = { ajo\So{Xj) , \bfo ) = So\fij)\bjo ) • (6.35) 

Under these assumptions the Darhoux matrix Ii5.16\) satisfies Ii6. 33\) and, 
moreover. 



aj$-i , \bj, ) = ^\bj) . (6.36) 



39 



Proof. We are going to verify that assumptions of the theorem imply fl6.33p . 
First, we will show that assumptions fl6.35l) imply (16.361) . Using (15.111) . (16.341) 
and (I6.35P we get: 

(6.37) 

IV ) = ^(/"/)IVo ) = $^(/i,)5'o(/i,)| Vo ) = ^h,) . 
Then, we compute 

{(ii'\bk') {aj\ hk)\jUk , / N 

Mfk' = , ' ^ = ^ " ^ , = -Xjf^kMjk , (6.38) 
Aj' — fik' A^/fc — -^j 

and, vice versa, Mjk = —Xji^k'Mjiy. Hence, 

1 1 
Kk'j' = r- Kkj , Kkj = r— Kk'j' . (6.39) 

Assuming J\f = I we proceed to compute ingredients of the formula (16.331) : 



^^0)^,_^M^l^ = j+^^JMMi^^ (6.40) 

j,k=i j,k=i 

where the second equality follows from: 



(primes can be dropped because we sum over the same set of indices). Then, 
we compute D{\~^) and decompose it into the sum of partial fractions: 



AT 



j,fc=l ^ j,k=l ^ 

Using (I6.40p . (16.360 and (16.390 . we get (after dropping primes): 

DiX-^) = DiO) + f (643) 
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Finally, 

~ 1 /X ^ ^\bi )Kiu( au\^^^ ^^^\bi )Wik( aul^^^ 
^D{X)<^-' = D{0)+J2 ' V \ + E X \ ^ ' 

(6.44) 

where 

AT TV 

Wjk = Kji{ ai \<l>^\bi )Kik = KjiMiiiXi - fii)Kik = {Xk - f^j)Kjk, 

i,l=l i,l=l 

(6.45) 

where we used $^ = J and (15.241) . Substituting Wjk into (16.441) and compar- 
ing the result with (I6.43p we get (I6.33p . □ 

Usually it is sufficient to assume Co(A) = = const (compare Re- 
mark l6.7l) and 5*0 = const (but the assumption So = I can be too restrictive). 

Proposition 6.15. We assume Sq = const, H = const, Sq = I, H'^ = H 
and SqHSq = H. We consider the Darhoux matrix \5.1b]) such that J\f = I, 
N = 2K and 

( «jo \ = {b]Q\H , ( aj+K,o I = ( &]+i^,o 1^, (6-46) 

( (^j+K,o \ = { Ojol'S'o , 'S'o|&j+A',0 ) = \ bjO ) , 

where j = 1, . . . , K , k = 1, . . . , 2K. Thus all these data can be expressed by 
{ aio !,...,( axo I ond Ai, . . . , Xk- The Darboux-Bdcklund transformation 
generated by such Darboux matrix preserves reductions: \l/^(A)i/\I'(A) = H 
and ^(A-i) = ^(0)^(A)^o. 

Proof: We apply Theorems 16.61 and I6.14[ It is enough to check whether the 
equations 

( «jo I = ( bio \H, { afo I = ( &]/o 1^' ( "i'o I = ( «io I'^o, ( &Jo I = ( 1*^0) 

are not contradictory. These equations imply ( ajQ\H~^ = ( ajo\SoH^^Sl. 
Hence, using Sq = I, we obtain the constraint S^HSo = H assuring the 
compatibility of both reductions. Finally, we denote j' = j + K. □ 
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7 Connections with other approaches 



In this section we shortly present some other methods of constructing the 
Darboux-Backlund transformation. We show how they are connected with 
the approach presented in this paper. 



7.1 Matrix- valued spectral parameter 

The name of Darboux first appeared in the context of the dressing transfor- 
mations in Matveev's papers (see for instance [15]) who extended the notion 
of Darboux covariance, known in the case of the Sturm-Liouville-Schrodinger 
spectral problems, on arbitrary differential operators [IB]. 

In order to apply Matveev's approach to Zakharov-Shabat spectral prob- 
lems (11. 2p the matrix spectral parameter is introduced: 

A := diag(Ai,...,A„) (7.1) 

(this notation should not be confused with the function A described in Re- 
mark 11.11) . We consider the linear problem of the form ^ |16] : 

= J] 5] Uuk.'^M^ + Kfc^A'^ , (7.2) 

j k=l k=0 

where Ui,kj and V^k are matrices which do not depend on Ai, . . . , A„ and 

1 1 



Mj := diag 



\l — CLj A^ — CLj 

The following theorem holds [6l |36] . 



Theorem 7.1. Equations are covariant with respect to the Darboux 

transformation 

^ = ^A-(T^, a = ^iAi^^\ (7.3) 



where is a fixed solution to (7.2) with A replaced by the diagonal matrix 
Ai = diag(Aii, . . . , A„i). 
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The linear problem fl7.2l) is closely related to the following special case of 
the standard Zakharov-Shabat linear problem fll.2p : 



j k=l ^ ^' fc=0 



Namely: 

vl/(A) = {$(Ai)pi,...,<l>(A„,K} , (7.5) 

where the notation used on the right-hand side (a matrix as a sequence of 
columns) is the same as in (12.131) and pi, . . . ,pn form a constant basis in C". 

The Darboux matrix generating the transformation ( I7.3P can be easily 
computed (using D = \E'\&~^). We get 

D{A) = $A^"i - ^lAi^^i . (7.6) 

If we put Ai = . . . = A„ = A, (i.e., A = A/), and pk = form the canonical 
basis in C" (i.e., {pi, . . . = J), then $(A) = ^(AJ) = ^(A). In this case 
we obtain 

D(A) = A/-^iAi^-\ (7.7) 

which is the starting point for the construction of the Darboux matrix by 
Gu and his collaborators [26l [271 [28l [811. Sometimes another form is used: 



D = J- A^iAr^^7\ (7.8) 
which is equivalent to (17.71) after changing A — > A^^. 

7.2 Transfer matrix form of the Darboux matrix 

A rational nxn matrix function -D(A), analytic at infinity, can be represented 
in the form [3], [25] : 

D{X)=Af + F{XlN-A)-'G , (7.9) 

where A is an x matrix. In is the unit matrix of order A^, and A/", F, G 
are matrices of sizes nxn, nxN and Nxn, respectively. Such representation 
is called a "realization" or a "transfer matrix representation" of D and the 
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number N (i.e., the order of A) is known as the "state space dimension" of 
the reahzation. Reahzations are not unique and can have different values 
of the number N. "Minimal realizations" have minimal value of N (and 
the minimal N is called the McMillan degree of D). Minimal realizations 
are unique up to a change of the basis in the state space (i.e., F FT~^, 
A TAT-^ and G TG, for some invertible x matrix T) [SI ES]- 



(7.10) 



Proposition 7.2. // ( [7.^ is a realization for D , then one of realizations for 
D^^ is given by 

D-\X) = AT-i - ^-^F{XIn -A + GJ\f-^F)-^GJ\f-^ 

The realization ^7.1(]\ ) is minimal iff ^7.9^ is minimal, see J2l 

The formula (17.101) can be verified by a simple but non-trivial computation. 
The obvious identity (A/tv - A + GM'^F) - {XIn - A) = GM'^F is very 
helpful, compare fl7.16p . 

Assuming JV = I we consider the so called transfer matrix 

Wa{x, A) = /„ - u;s~\A - xi^r'n, , (7.11) 

where A, S", 111,112 are some matrices (the star denotes a matrix conjugate, 
but this is not very important at this moment) and, moreover, the following 
operator identity holds: 

AS-SB = UiU;. (7.12) 

Matrices A, B,]li,]l2, S satisfying (17. 121) are said to form an S-colligation 



The transfer matrix (17.111) can be used to generate solutions to integrable 
systems by the Darboux-Backlund transformation, see [621 163] . We can make 
the following identification: 

S = M , s-^ = k. 



A = diag(Ai, A2, . . . , Xn) , B = diag(/ii, /i2, • • • , IJ-n) 
and, finally 

( Cio I 

u; = {\h),\b2),...,\bN)) , ni= \ ' 

\ ( CtAT I J 



(7.13) 



(7.14) 
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Corollary 7.3. The symmetric representation of the Darhoux matrix Ii5.16\) 
can be identified with the transfer matrix form ( [7. where A is diagonal. 
The identity ( [7. coincide with the definition Ii5.13\) of the matrix M . 

Constant matrices A of more general form correspond to generalizations 
of (15. ip (multiple poles are allowed). 

In order to show a flavour of the transfer matrix technique we present one 
of typical results. Note that the proof of Proposition 17.41 is similar to some 
steps in the proof of Theorem 15.51 



Proposition 7.4. We assume the identity jl.l^ and D = wa is given by 
the formula (TJl^. Then 



D-' = In + U;{B-\lN)-'S-'ni (7.15) 



Proof: We compute 

(/„ - u;s-\A - A/jv)-^ni) (/„ + u;{B - xir,)-'s-'u,) = 4+n;xni , 

where X is x matrix given by 

X = {B-X)-^S-^-S-\A-X)-^-S-\A-X)-\AS-SB){B-X)-^S-^ . 

Now, using the obvious identity 

AS - SB = {A- X)S - S{B - X) , (7.16) 

we decompose the last component of X into the sum of two terms which 
immediately cancel with the first two components of X. Therefore X = 
which ends the proof. □ 

Vectorial Darboux transformations constitiue one more approach to Dar- 
boux transformations, applied mostly in 2 + 1-dimensional case [151 Bl] . Al- 
though this technique needs no analogue of the Darboux matrix but the 
Darboux transformation is expressed by a Cauchy-like matrix and important 
role is played by operator identities like (I7.12p . Comparing the results of [H] 
and [62] we conclude that both methods are in a very close correspondence 
(note that the matrix S of [52] corresponds to the matrix $ of [H]). 
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7.3 Factorization in loop groups 

Given a Lie group G we define the loop group of G as tlie group of smootli 
functions 7 : 5*^ — ^ G, wliere denotes tlie unit circle on the complex plane 
(|A| = 1) [29l [59]. An important role in the loop group theory plays the 
Birkhoff factorization theorem. The Birkhoff decomposition is closely related 
to the Riemann-Hilbert problem which provides a rigorous background for 
the inverse scattering method [56], see also [29] . 

In general the Birkhoff factorization is not explicit. The explicit cases 
are closely related to the construction of Darboux matrices [751 HSl CZj (and 
also to the construction of finite gap solutions), compare similar ideas in the 
soliton theory [33 ESj . The approach based on the so called cc-ideals is one 
more link between the loop group theory and the theory of solitons [531 El] • 

From geometrical point of view the Lax pair consists of commuting differ- 
ential operators and their compatibility can be interpreted as the condition 
that a one-parameter family of connections is flat: 

[9i-f/i(a;,A), 92-f/2(x,A)] = (7.17) 

(f/i, U2 are matrices depending on x through some fields, say u). The "trivi- 
alization" E oi a. solution u is defined as the solution of the system: 

E,,= -EU,, E{0,\) = I. (7.18) 

Then E{x,X) is holomorphic for A G C, see [THj. The function E{x,X) is 
also referred to as an "extended solution" , an "extended frame" or simply a 
"frame". Comparing fl7.18p with (11.21) we can identify E = Actually, 
fl7A8|l is the adjoint of ([L2D, see also 

Theorem 7.5 (Birkhoff). The multiplication map fi 

fi : L+{GL{n,C)) x L_(GL(n,C)) ^ L(GL(n,C)) 

is a diffeomorphism onto an open dense subset of L{GL{n,C)) , where 

• L^{GL{n, C)) is the group of holomorphic maps : C ^ GL{n, C) 

• L„{GL{n,C)) is the group of holomorphic maps h_ : Ooo —>■ GL{n,C) 
such that /i_(cxd) = /, where Ooo is a neighbourhood o/A = 00. 

• L{GL{n, C)) is the group of holomorphic maps from OoqCiC to GL{n, C) 
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Corollary 7.6. Suppose that h_h^ lies in the image of jj,. Then, by virtue 
of the Birkhoff theorem, there exists a unique pair f± G L±{GL{n,C)) such 
that h^h+ = f+f-. One can interpret it as a "dressing action" of on hj^ 
and /+ is the result of this action, which is denoted by /?._tl/i+ = /+. 

The dressing action seems to "forget" about /_. However, it is worthwhile 
to stress that this is /_ which should be identified with our Darboux matrix. 
On the other hand the element /i„ is deeply hidden (almost non-existing) 
in other approaches to the construction of Darboux matrices. In order to 
explain the dressing action generated by the Birkhoff decomposition we will 
present the binary Darboux transformation (12.11) in the framework of the 
loop group approach, following [76j. 

We assume that E{x,X) G L_|_(GL(r;,, C)) is given, and we choose the so 
called "simple element" h\^^^^ .„ G L_{GL{n,C)): 

hx,,,U>^) = I + ^^^n , (7.19) 

where Xi, fii are complex parameters, and vr is a constant (x-independent) 
projector in C" (i.e., vr^ = vr). One can easily see that h'^^^^^ = /i^i,Ai,7r, 
compare ([2A]) and (Q- 

Then, the Birkhoff theorem states that there exists E G L^{GL{n,C)) 
and D G L_(GL(n, C)) such that 

hx„f,„nE{x, A) = E{x, X)D{x, A) , (7.20) 

provided that the product on the left-hand side belongs to certain "open 
dense set" of L{GL{n,C)). Now, both the exact form of D and this "open 
dense set" can be found by direct calculation. It is sufficient (similarly as in 
all other approaches discussed earlier) to compare the residua at both sides 
of the equation (I7.20p . Hence 

D = I + ^1^P, (7.21) 
A — Ax 

where P is defined by (12.111) . where Vker = kervr and Vim = im vr. The open 
dense set from the Birkhoff theorem is defined by: kerP fl imP = {0}. We 
remark, by the way, that the Birkhoff theorem assumes the isospectral case 
and the canonical normalization {J\f = I) . 

Note that (17:201) implies E = hED~^ = h'^-^D-^ = (D^h-^)-^ (where 
h denotes the simple element). Therefore, = E^^ = D'^h^^, which is 
equivalent (because h does not depend on x) to the usual formula ^ = D"^ 
(compare Remark 11.41) . 
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8 Invariants of the Darboux transformation 



The Darboux transformation changes the matrices U,y into new matrices Uk 
of the same form. By invariants of the Darboux transformation we mean 
constraints on coefficients of f/^ which are preserved by the transformation, 
see [11] (compare also [69j, where one may find many examples). The in- 
variants are very useful in the construction of Darboux matrices in purely 
algebraic way, without referring to the special boudary conditions and to the 
scattering theory (which is a usual practice, compare [271 [76t [8l]). 

Here we simplify the approach of [H] and extend it on the non-isospectral 
polynomial case. Moreover, we show that our approach works also in much 
more general case: when the Lax pair is singular at some fixed values of the 
spectral parameter. In this section we denote Ui = U, U2 = V. 



8.1 Linear invariants for polynomial Lax pairs 

We consider Lax pairs with the following A-dependence: 
00 00 
U = Y, UkX""-" = X''u, V = J2 = ^^'v > (8-1) 

fc=0 fc=0 

where N, M are fixed positive integers (not to be confused with the notation 
of previous sections) and Uk = Uk{x), Vk = Vk{x) (A; = 0, 1, 2, . . .). Usually the 
sums are finite (i.e., Uk = Vk = for sufficiently large k), and this typical case 
(polynomial in A and A~^) corresponds to many classical soliton equations. 
In particular both U and V can be polynomials in A (in this case, for N = 1, 
we get famous AKNS hierarchy). We also assume a similar A-dependence of 
the derivatives of A: 

00 00 
A,i = , \2=Y. bkX'''-' , (8.2) 

where N', M' are given integers fixed by the assumption ag 7^ 0, 60 7^ (in 
the nonisospectral case). The coefficients ak = ak{x), bk = bk{x) have to 
satisfy compatibility conditions resulting from A, 12= A, 21 (some examples 
can be found in [IH [H EH)- 

We consider the Darboux transformation of U + HV, where H = H{x, A) 
is a fixed function 

H{x, A) = A^-*^/i(a;, A) = X^-^'{ho + hX-' + h2X-^ + ...), (8.3) 
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where Hq, hi,h2, . . . are given functions of x. We assume that H is unchanged 
by the Darboux transformation (and U, V are transformed, as usual, accord- 
ing to (11.41) ). The Darboux transformation yields 

{ij + HV)D = Da+hD,2+D{U + hV) , (8.4) 

which reduces to 

{u + hv)D - D{u + hv) = \-^D,i +h\-^^Da ■ (8.5) 
We assume that D is analytic at A = oo : 

D = To + Ti\-^ + + . . . , det To ^ , (8.6) 

i.e., D = X^'^D, where D is given by (13.11) . 

The idea of linear invariants is quite obvious. Suppose that for A ~ cxo 
the right-hand side of (18.51) behaves as A~^, where K ^ 1. Then the first 
K terms of the Taylor expansion (in A~^) of the left-hand side are equal to 
zero. The first two of these equations read: 

(uo + hovo)To = To{uo + hovo) , 

(ui + hovi + hiVo)To + (uo + hoVo)Ti = To{ui + hovi + hivo) + ri(no + hoVo) . 

The asumption uo+hoVo = implies ■Uo+^o'^o = (provided that det To 7^ 0). 
Then, adding the second assumption: Ui + h^vi + hiv^ = 0, we obtain as a 
consequence ui + h^vi + hiv^ = 0. Thus we have two expressions invariant 
with respect to the Darboux transformation. Considering the first k (where 
k ^ K) equations we get an invariant system of k equations. 

We proceed to estimate K. The leading terms of the right-hand side of 
(18.51) are given by: 

A-^(To,i -aoTiA^'-^ + ...) + ho\-^'{To,2 -hTA''''^ + . . .) (8.7) 
Therefore K > kmaxi, where 

kmaxi = -1 + min{iV, M, N + 2 - N' , M + 2 - M'} , (8.8) 
what can be summarized as follows. 
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Proposition 8.1. Suppose that ^ k ^ kmaxi o,nd Hq, hi, . . . ,hk are given 
functions of x. Then the system of k + 1 linear constraints 



+ 5^M,_, = 0, (j = 0,l,...,A;) , (8.9) 

i=0 

is invariant with respect to Darboux transformations such that detTo ^ 0. 

In some special cases, we can formulate stronger propositions (i.e., we have 
more invariants). In the isospectral case we can replace kmaxi by 

Ca.i = -l + min{iV,M} , (8.10) 

(the same result is valid when N' ^ 2 and M' ^2). In the case of the 
canonical normalization (Tq = /) we can replace k^axi by 

Caxi = min{A^, M,N + 1- N',M + 1- M'} . (8.11) 

Below we present one more example. 

Proposition 8.2. Suppose that min{M, N + 2-N',M + 2- M'} > N 
(it implies, in particular, kmaxi = N — 1), functions Hq, hi, . . . , hk (k ^ N) 
are given, and Tq assume values in some matrix Lie group G. Then, the 
following system of k + 1 linear constraints is invariant with respect to the 
Darboux transformation: 



+ Y,hiVj~i = 0, {j = 0,l,...,k-l] 



(8.12) 

«fc + ^ hiVk-i e g , 

1=0 

where g is the Lie algebra of the Lie group G. 

The proof of this proposition is analogical to the proof of Proposition 18. 11 
we consider coefficients by powers of A"^ in (18. 5p . Only the last step has to 
be treated in a different way. Assuming that the first k constraints hold, the 
coefficients by A~'^ yield 

k Z' ^ \ 

Uk + ^ hiVk-i = To I Mfc + ^ hiVk-i j Tq^ + (5fcjvTo,i T^^ . (8.13) 

Now the proof follows immediatelly from well known properties of matrix Lie 
groups (TgT~^ C g and T,i G g, provided that T = T(x) G G). 
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8.2 Bilinear invariants for polynomial Lax pairs 

Assuming the polynomial form (18.11) of U, V we consider the Darboux trans- 
forms of bilinear forms Tr(f/^), Tt{V^) and Tt{UV). We present computa- 
tions for the last case (the other two cases are analogical). In this setion we 
use notation: A ■ B = Tt{AB). From (ll.4p we get 

Tr{UV)-TT{UV) =TT{DaD-^D,2D-^ + D,iVD-^ + D,2UD-^) . (8.14) 

The leading terms of the right-hand side of (I8.14p read 

^-(7v+M)Tr ((To,! -aoT^X^'-' + . . .)T^\T,,2 -boT,\^'''' + . . .)T,-') , 

A-^Tr ((To,! -aoTiA^'-2 + . . .)voT,') , 

A-^^Tr ((To,2 -6oTiA*^'-2 + . . .)uoT,~') . 

(8.15) 

Thus the right-hand side of (I8.14p behaves as A""^, where K will be estimated 
below. 

We assume that D and D^^ are analytical at A = cxd (i.e., detTo 7^ 0). 
Considering coefficients by A^-' (j = 0,1,2,...) in the formula (18.141) . we 
obtain the following invariants: 

/o := Uo-vo , 

fi := Uo ■ Vi + Ui ■ Vo , 

f2-=Uo-V2 + Ui-Vi + U2-Vo, (8.16) 



fk := Uo-Vk + Ui- Vk-1 + . . . + Uk ■ vo , 
where k < K. In order to formulate a more precise statement, we define 

where = M + iV - 1 + min{0, 2 - N' ,2 - M' ,4 - M' - N'} and kmaxi is 
given by (18. 8p . 

Proposition 8.3. Bilinear expressions fk (k = 0, . . . , kmax2), given by l{8.16\) . 

are preserved by the Darboux transformation (i.e., fk = fk) provided that 
detTo ^ 0. 



51 



Remark 8.4. If M > N ^ 0, N' ^ N + 2, M' ^ M + 2, then kmax2 = k^axi- 

In some cases we can formulate stronger propositions. For N' ^2, M' ^ 2 
(including the isospectral case) kmax2 in Proposition 18.31 can be replaced by 

k'max2 = -1 + min{iV, M, N + M] . (8.18) 

If the normalization is canonical (Tq = /) we can replace kraax2 by 

kmax2 — ^i^{^maxl5 ^mn} 5 (8.19) 

where k'^^ = M + N + min{l, 2 - A^', 2 - M', 3 - M' - A^'}. 

Analogical considerations can be done for Trt/^ and TrV^. To obtain the 
final results (see below) it is enough to substitue M —>■ N, M' N' in the 
first case, and A^ — > M, A^' — > M' in the second case. 

Proposition 8.5. Suppose that ^ k ^ kmaxs, where 

k^axs = mm{N - 1, N + 1 - N' ,2N - 1,2N + 1 - N' ,2N + 3 - 2N'} 
and go, gi, . . . , g^ are given functions of x. Then the bilinear constraints 

go := uo-uo , 

gi := Uo ■ ui + Ui ■ Uq , 

g2 := Uo ■ U2 + Ui ■ Ui + U2 ■ Uq , (8.20) 



9k ■= Uo-Uk + Ui- Uk-i + ... + Uk-Uo , 
are preserved by the Darboux transformation such that detTo ^ 0. 

Proposition 8.6. Suppose that ^ k ^ kmaxi, where 

kmaxi = mm{M - 1, M + 1 - M', 2M - 1, 2M + 1 - M', 2M + 3 - 2M'} 
and ho, hi, . . . ,hk are given functions of x. Then the bilinear constraints 

ho ■■= vo-vo , 

hi := t;o ■ t^i + t^i ■ t^o , 

/i2 := f ■ + 'Wi ■ 'Wi + ^"2 ■ f , (8.21) 



hk ■■= Vo-Vk + Vi- Vk-i + ... + Vk-Vo , 
are preserved by the Darboux transformation such that detTo 7^ 0. 
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8.3 Invariants for general Lax pairs 

Let us consider matrices U and V in the neighbourhood of A = Aq, where 
U, V have poles of A^-th and M-th order, respectively, i.e., 

oo oo 

f/ = ^Mfe(A-Ao)'-^, V = Y,Vk{X-Xo)'-^' , (8.22) 

k=0 k=0 

where Uk = Uk{x), vu = Vk{x) {k = 0, 1, 2, . . .). We are going to show that 
the general case reduces to the polynomial case discussed above. Indeed, it 
is sufficient to use the following parameter z in the neighbourhood of Aq: 

z = {\-\o)-\ (8.23) 

and then the Lax pair (18.221) becomes identical with (18.11) . Note that z ^ oo 
for A — > Aq. We assume that the Darboux matrix D is analytic at Aq: 

D = To + Ti(A - Ao) + T2(A - Ao)' + . . . = Tq + z'^T, + z'^T^ + ... 

where matrices depend on x. In the nonisospectral case we transform the 
equations (II. 3p to the form (18.21) : 

z,,= -z^L,{x,Xo + z-^) , (8.24) 

where L^, have to be expanded in the Laurent (or Taylor) series at z = oo. 

In order to obtain linear invariants we consider the linear combination of 
matrices U, V, given by: 

U+{X-Xo)^'-''hV (8.25) 

where 

oo oo 

h = h{x, y- A) = ^(A - \^fhu{x, y) = Y, ^kZ-" (8.26) 

k=0 k=0 

is a given scalar function, holomorphic at A = Aq. Finally, we arrive at an 
exact analogue of Proposition 18. 1[ 

Bilinear invariants can be treated in the same way. We obtain exact 
analogues of Propositions 18.31 18.51 and 18.61 

Corollary 8.7. The polynomial case can be treated as a special subcase, 
defined by Xq = oo. It is enough to change variables in formulas (E! 
A — >• A~^ (and Xq Xq^). Then, making the limit Aq — > 0, we get l\8.1]) . 
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8.4 Application to the KdV equation 

We will show advantages of Darboux invariants considering the case of the 
KdV equation. Our approach consists in characterizing the Lax pair in terms 
of some algebraic constraints (see fiSi [Ti] ) and then showing that these con- 
straints are preserved by the Darboux-Backlund transformation. 

Proposition 8.8. The Lax pair ( [^. can be uniquely characterized by the 
following set of alebraic constraints: 

1. U is linear in X (U = UqX + Ui), Tr f/ = , 

2. V is quadratic in X (V = VqX'^ + ViX + V2), TiV = , 

3. uq = [ ^c\]' ui is off-diagonal. 



4- Uq — = , Ml — G g , where g is the 1-dimensional Lie algebra 
generated (y^^ q 

5. Uq ■ Vi + Ui ■ Vq = , Ml ■ Ml + 2mo ■ M2 = , 



6. U{X) = U{X), V{X) = V{X) . 
Proof: The first four properties imply the following form of U, V: 
\ f P 



^° = t -1 j ' - V « 

00\ f Ap\ f -a h 



= > -4 j ' = U j ' \ c a 



(8.27) 



where M,p, g, a, 6, c are some complex fields. Bilinear constraints O yield 

-8p = -8 , 8pg - 86 = , (8.28) 

i.e., p = 1, q = b. Now compatibility conditions yield the KdV equation 
f l4.10p and expressions fl4.14p for a, b, c. The last property implies m G M. □ 

The first two constraints are preserved by any Darboux transformation 
constructed in the standard way, e.g., using Corollary 13.71 (and the traceless- 
ness is preserved by virue of Remark 13.31 provided that det A/" = const). The 
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constraints E] impose restrictions on and p^, see Section 16. 2[ In order to 
preserve the third constraint we have to use freedom in the choice of the nor- 
mahzation matrix Tq = M . From the first equation of fl3.19p we get (taking 
into account the form of Mq given by the third constraint): 

Af=/(^;). (8.29) 

where /, a are some functions. In the sequel we put / = 1 (thus det A/" = 1). 
Then, denoting Ti = ( '"^ '"^ ] , we rewrite the second equation of (13.191) 



as 



C3 C4 

" M = f -\ M + f " « M.(8.30) 

u u J y u — a a J \ ci — ac2 — C4 C2 y \ a,i y 

Hence: 

a = — C2 , u = u — a"^ — ac2 + a,i +Ci — C4 . (8.31) 

The constraint m is preserved by virtue of Proposition 18. 2[ Other propo- 
sitions from Section [8] are too weak for our present purposes. Indeed, in the 
KdV case we have kmax2 = and kmaxi = 1- Therefore the preservation of 
the constraints does not follow from Propositions 18.31 and 18. 6[ 

Fortunatelly, the special form of matrices uo,vo and Tq (given by (18.291) 
with / = 1) for KdV equation allows us to reconsider the behaviour of leading 
terms (18.151) . We easily see that any matrix product containing only matrices 
from the set {Tq, Tq"^, Tq,i , Tq,2 , Uq, Vq} (and among them at least one matrix 
from the set {To,i , To,2 , Uq, Vq}) is proportional to Uq, and, as a consequence, 
it has vanishing trace. Hence, in the case of the KdV equation the Dar- 
boux transformation preserves constraints (18.161) for k = 0,1 and constraints 
(I8.2ip for k = 0, 1,2. In particular, the constraints [5] of Proposition 18.81 are 
preserved. 

Corollary 8.9. The Darhoux transformation (defined as in Corollary 3.1) 



preserves all constraints defining the KdV Lax pair (see Proposition \8.t^) 
provided that we impose reality restrictions on Xk, Pk (see Section W^) and 
fix the normalization matrix according to the formula h8.2^^) where / = 1 and 
a is expressed by the matrix Ti, namely a = —02- 

In the case of the elementary Darboux matrix det Tq = and considera- 
tions presented in this section are not apphcable. It would be interesting to 
extend the theory of Darboux invariants on the case detTg = 0. 
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9 Concluding remarks 



In this paper we gave a unified view on the Darboux-Backlund transforma- 
tions for 1 + 1-dimensional integrable systems of nonhnear partial differential 
equations. In particular, we discussed in detail relationships between various 
approaches to the construction of the Darboux matrix. 

Darboux-Backlund transformations have been extended in many direc- 
tions. First of all, they are applicable to 2 -|- 1-dimensional integrable systems 
P [21 |27l Sg EH , including self-dual Yang-Mills equations [271 [551 [78] . Then, 
we have 0-|- 1-dimensional systems, e.g., ordinary differential equations of non- 
linear quantum mechanics [211 [211 [SZ] . Darboux transformations were also 
constructed in the supersymmetric case [121 [IB] and in the non-commutative 
case [65,]. 

Matrix representations of spectral problems and Darboux transformations 
are not always convenient. Impressive examples are associated with Clifford 
algebras. It is enough to compare the paper [TS], where mainly the matrix 
approach was used, with subsequent papers [SI HZ], which are much shorter, 
more general and more elegant. All these papers consider binary Darboux 
transformation. An extension on multipole case is not so obvious, compare 
|20j . where some progress in this direction is described. There exist other 
generalizations of the Darboux transformation on spectral problems with 
values in abstract associative algebras [Tni[21|- 

The discrete case is (to some extent) very similar to the continuous case. 
Many aspects (e.g., those concerning the rational dependence on A and the 
loop group structure) are just repetitions from the continuous case, compare 
[211 [Ml [ISj- It is tempting to apply the ideas of time scales [S], all the more 
so that in the "classical" case of the pseudospherical surfaces we succeded to 
construct the Darboux-Backlund transformation on arbitrary time scale [19] , 
thus treating the discrete and continuous case in a uniform way. However, 
some points seem to be more difficult in the discrete case, e.g., Darboux 
invariants are not formulated yet. Actually, it is not so easy even to find 
an appropriate discretization of a given integrable system, especially if the 
associated linear problem is non-isospectral. 

Acknowledgements: I am grateful to Maciej Nieszporski for many fruiful 
discussions. 
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